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Abstract 

We extend an IR-deformed soft- wall AdS/QCD model at zero temperature to a model at finite 
temperature and perform hydrodynamics. To have the in-falling boundary condition to make the 
hydrodynamic analysis possible, we treat the infrared energy scale factor in our metric as a temperature- 
depending parameter. Then by carrying out the hydrodynamic analysis, we compute the transport 
coefficients, the diffusion constant and the shear viscosity, through the linear response theory. 



1 Introduction 



To carry out the analyses in the strongly coupled region would be regarded as the key of many important 
unsolved problems in the contemporary high energy theoretical physics. One of these problems will be the 
low energy dynamics of quantum chromodynamics(QCD). Actually, many ideas in the non-perturbativc 
analysis of QCD have been motivated. Among these, there are three interesting classes, chiral dynamic 
model based on large N c [1], the lattice gauge theory [2] and holography [3,4] (gauge/gravity correspon- 
dence, AdS/QCD, etc). 

In the lattice gauge theory, Minkowski space-time in the original theory is replaced with a finite volume 
Euclidian discrctizcd space-time with analytic continuation of time direction to imaginary time direction. 
As a result, the degree of freedom of theories becomes finite, and the non-perturbativc numerical analysis for 
the action itself, which is Monte Carlo simulation, becomes available. However, Monte Carlo simulation 
has a notorious problem named as the sign problem, when fermions are involved in the Monte Carlo. 
Currently, no mean for it has been invented yet, and the current analyses in the lattice are always carried 
out by getting around it [5]. 

On the other hand, it is well known that the holography is a duality between the strongly coupled field 
theories and the weakly coupled gravities. It actually originates in the superstring theory in which the 
quantized open and closed strings at low energy can be identified with particles known in field theories 
and gravities [6]. As a result, behaviors of the low energy open strings describe the supersymmetric gauge 
theory. Then as an important matter in the correspondence, U(N) supersymmetric gauge theory arises 
on N overlapped D-branes from the open strings sticking to these N overlapped D-branes at low energy. 
On the other hand, D-branes can be identified with black branes in the supcrgravity [7] . 

Thus one has two ways to describe a low energy D-branes. As a result, it is known that one can 
conjecture the duality between p-dimensional large- N supersymmetric SU (N) gauge theory with large 't 
Hooft coupling and near- horizon geometry of black p+ 1-branes in the condition that the quantum effect of 
gravity and the length of string can be neglected. Although there is no exact proof for this correspondence 
at this moment, it is particularly expected that the duality between the M = 4 four-dimensional large- TV 
SU(N) supersymmetric gauge theory [3,4] and the five-dimensional anti de Sitter space is valid. 

One of the great advantages in the gauge/gravity correspondence compared with the lattice gauge 
theory would be that it is irrelevant to the problem arising when one involves fermions like the sign 
problem. Because the space-time in the gauge/gravity correspondence is not discretized space, and the 
main analysis are carried out in the weakly coupled gravity side. 
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However, the current gauge/gravity correspondence also has a problem. It is that the dual gauge theo- 
ries are always non-realistic as long as the gravity side is a solution (Top-down model), and gauge theories 
become supersymmetric. Due to this, many results in the field theory side in the current gauge/gravity 
correspondence is the ones independent of detail of theories or no more than qualitative ones or results just 
in supersymmetric model. On the other hand, once getting away from the study based on a solution in 
the gravity side, constructing holographic models in bottom-up way is also conducted energetically [8-10]. 

Anyway, the point that the gauge/gravity correspondence can be irrelevant to the peculiar problem 
in the treatment of fermions in the lattice gauge theory would be one of the great advantages. For this 
reason, the low energy dynamics of QCD has recently been studied intensively in the framework of the 
gauge/gravity correspondence. This came to the motivation on an IR-deformed soft-wall AdS/QCD model 
□ [11], which is a kind of the holographic bottom- up model [9]. 

Recently, we have made the extension of an IR-deformed AdS/QCD model [11] to the finite temperature 
system [12, 13], where the deformed bulk vacuum and potential term have to be introduced for the scalar 
field to satisfy the equations of motion. This is because if one straightforwardly extends the model to the 
finite temperature, the solution of the equation of motion diverges for the dilaton. Only by deforming 
the bulk vacuum and potential term, one can obtain the smooth solution for the dilation. With such 
a treatment, we have examined the critical temperature of chiral symmetry breaking [12] through the 
analysis of the quark number susceptibility, and the meson spectrum [13]. In the analysis of the meson 
spectrum, we have carried out the numerical analysis of equation of motion for the fluctuations on the 
bulk gravity. Such a numerical analysis for the mass spectrum has been considered as a basic method of 
analyses in holographic QCD as well as the hydrodynamics. 

One of the crucial points of the hydrodynamics is that it can be considered to be independent of 
detail of theories, and so independent of the transport coefficients, where the transport coefficients are 
the ones defined in the framework of the hydrodynamic analysis. In particular, the ratio rj/s between the 
shear viscosity (77) and the entropy density (s) is the quantity characterizing the actual QCD. For this 
reason, the ratio rj/s has been examined very much in the gauge/gravity correspondence [14]. Besides, 
the holographic hydrodynamics has played an important role in the longstanding problem in the causal 
hydrodynamics [15]. and quark-gluon plasma described by Bjorkcn flow [16]. 

With such a circumstances, turning to the hydrodynamics, we are going to work out in this paper the 
holographic hydrodynamics in an IR-deformed AdS/QCD model at finite temperature studied in [12,13]. 
Some interesting features mentioned above will be investigated in the future. 

*We have referred to this IR-deformed soft-wall AdS/QCD model as a predicted model in our former papers [12, 13]. 
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Then wc would like to mention the organization of this paper. In ChapJ^l we are going to review the 
hydrodynamics and the transport coefficients obtained from linear response theory. In ChapJ3] we will 
introduce our holographic model, and show the way how the model is extended to the finite temperature 
system. It is shown that if one takes the same way as the one in ref. [13], it turns out that the analysis 
becomes too complicated to be carried out. Therefore, in this paper, we will propose an alternative way, 
which is simpler than the one in ref. [13] as the quark-gluon plasma is considered to be above the critical 
temperature. It will be seen that there are in general four options in the numerical calculations for some 
factors, while only one of them is physically acceptable. In Chapf?] we will sort out the notation used in 
this paper. In Chap|5] we will carry out the hydrodynamic analysis for the fluctuation of the balk gravity 
dual to U(l) baryon current in the scalar mode. In the analysis, a necessity to cancel the divergence at the 
horizon arises as usual. As our model involves the dilaton which makes the analysis complicated. However, 
we will show how to contain it by exploiting the integral constant. By using GKP-W relation [4], we will 
read out the retarded Green function for U(l) baryon current in the scalar mode and its diffusion constant. 
In Chap|51 we are going to carry out the hydrodynamic analysis for the fluctuation of the balk gravity dual 
to SU(2) flavor current in the vector mode as well as in Chap|5] Again by using the GKP-W relation [4], 
we will read out the retarded Green function. In ChapfT] we shall carry out the hydrodynamic analysis 
for the fluctuation of the tensor mode in the balk gravity dual to the energy-momentum tensor. Thus by 
using the GKP-W relation [4] , we can directly read out the viscosity with the retarded Green function and 
evaluate the ratio rj/s. Our summary and conclusions will be presented in Chap[5] 

2 Brief Review on Hydrodynamics and Linear Response Theory 

To begin with, we would like to describe the basic matters in hydrodynamics and linear response theory 
used in this paper. The description in this chapter is basically following the review papers [17]. 

The hydrodynamics is an effective theory to describe the macroscopic dynamics at large distances and 
time-scales. Conserved quantities are considered to survive in such large distances and time-scales, and 
the energy-momentum tensor T^ u is one of conserved quantities. The hydrodynamics is formulated by the 
hydrodynamic equation instead of action principle. The hydrodynamic equation for the energy-momentum 
tensor is given as 



d^ v = 0, 



(2.1) 



where 



= (e + P)u^u v + Pg^ + r 



(2.2) 
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with energy density e, pressure P, local fluid velocity w' 1 and the given as 



t"" = -n ( d"u v + d v u" - \^ v d a u a ) - (V^d a u a . (2.3) 
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Here r\ and £ mean shear and bulk viscosities, respectively. In a curved space, it is given as 



(2.4) 



where P MI/ = + 1*^". In the above, 77 and £ are regarded as the transport coefficients. In what follows 
we consider the fluid in the rest frame, = (1, 0, 0, 0). 

In eq. (|2.2[) . substituting u M = (1, 0, 0, 0) and expanding as = r\^ iv + h^ v (rj^ — diag(— 1, 1, 1, 1) and 
h^ u means fluctuations), it turns out that 8{r xy ) is given as 

S(T xy ) = iujrih xy , (2.5) 

where we have represented the formula in the momentum space. It is known from the linear response 
theory that the response of an operator O for the external field <f>^(k) is given as 

S(0(k)) = -Gg°0(k), (2.6) 

/oo 
d 4 xe~ lkx ([0{t, x), O(0, 0)])0(t) is the retarded Green function. Then by compar- 
-00 

ing eg. (12.51) with ca. (|2.6[) . we can obtain the Kubo formula with regard to the shear viscosity as 

77 = - lim -ImG", (2.7) 
' mow y ' 

/OO 
d 4 xe- tkx ([T x y(t, x), T xy (0, 0)])0(t). 
-00 

Next we turn to the diffusion in the conserved current and the energy-momentum tensor. Toward a 
conserved current j M with = d^j 11 , the constitutive equation for the conserved current is given as 

j = P u- DVf = -DVj\ (2.8) 

where j = p and D mean the charge density and the diffusion constant respectively, and we have taken 
account of the fluid rest frame. Then the Fick's low is valid as 

dtp - DV 2 p = 0. (2.9) 

This gives the following dissipation relation for the charge density, 

uj = -iDk 2 , (2.10) 

/OO 
d 4 xe~ tkx 
-00 

([p(t, x), p(0, 0)])6(t). Next, we consider the diffusion in the the energy-momentum tensor T tl with i = x, y. 



First, from the constitutive equation in the curved space, one can obtain 

T zl = -—l^d z T u = -i]d z u\ (2.11) 

Here, we have arranged the k along with the x 3 -axis as k = (0, 0, k) for simplicity. Then, from d^T^ 1 = 0, 
one can obtain 

d t T u ^—d 2 T u = 0. (2.12) 

e + P z y J 

This gives the dissipation relation in the energy- momentum tensor T tl , 

u = -i—!]—k 2 , (2.13) 
e + P ' v ; 

which will be the pole in the retarded Green function u . 

Finally, we write down the decomposition of the energy-momentum current and the U(l) current under 
the little group 5*0(2) toward i-direction (i = x,y): 

scalar mode: T 00 , T 03 , T 33 , T\ and J , J 3 , 
vector mode: Tq^, and Ji, 
tensor mode: — S^T \/2. 

3 An IR-deformed AdS/QCD Model at Finite Temperature 

We will start with the following geometry, which is deformed from Schwarzschild AdS§ black hole geometry 
in the IR-region by a factor /j, g as 

ds 2 = a 2 (z) {-f{z)de + £ dx* + (3.1) 

with a 2 (z) = (1/z 2 + H 2 g)/l 2 and f(z) = 1 — (z/zo) 4 (zq means the location of the horizon) and we have 
put the AdS radius I as 1 in what follows. The coordinate z is in the relation with the usual coordinate r 
as z = 1/r. This geometry has Hawking temperature T = 1/(jtzq), and is asymptotically AdS§ space-time. 



We will consider the following U(l) x SUl{2) x SUr(2) soft-wall model with scalar field on the back- 
ground (|3.ip as 



I d 5 x V^ge-^ ( R + ^ ) + I d 5 x V^ge^ 



1 

~2 fMN-F 



4 5u ( i 



(i) 



+Tr 



^—(Fl.mnFI in + F R , MN F% N ) + \D M X{z)\ 2 m 2 x \X(z)\ 2 - ±\X(z)\ 4 

SU(2) 4 



,(3.2) 
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where m\ = —3 [8,9] and / = 1. The trace is performed for the SU(2) algebra mentioned in what follows. 
We write Fa/at, F l>mn and F R>MN as F MN = d^A y -dyA^ and F l<M n = d^B LtI/ -d v B Lill -i[B LtI/ ,B Ltll ], 
where B^ ^ = £?£ ^t 0, with the SU{2) Lie algebra t a (a — 1, 2, 3), and now we have skipped to describe for 
R part. As for the dual operators for these in the boundary theory, for example see the table in ref. [8]. 
Using these, the covariant derivative can be written as DmX = 8mX + i(B^ .mX — XBr m)- 
The bulk vacuum of the scalar field is given as 

X(z) = ^ 1 2 , (3.3) 

where I2 means 2x2 unit matrix, and v(z) is given in TablcfTJ As it is known that v(z) behaves around 
the boundary as 

v(z)=m q (z + ^z 3 + 0(z 5 ). (3.4) 

From AdS/CFT correspondence, m q and a can be interpreted as the quark mass and quark condensate, 
respectively. As for £, see Tabled] and [2] 



Model 


v(z) 


Parameters 


lib 


z(A + Bz 2 )(l + Cz 4 )- 5 ^ 8 


A = m q (, 


B^a/C C = (£ 2 /(/^7 2 )) 4/5 



Table 1: These are taken from ref. [11]. The numerical values for these are given in Table(2j which are fixed 
by minimizing the breaking of the Gell-Mann-Oakes-Renner relation, /| m% = 2m q a, at the 1% level and the 
experimental values: = 139.6 MeV and f„ = 92.4 MeV. 



Model 


A 


m q (MeV) 


as (MeV) 


7 


(J-g 




lib 





4.07 


272 


0.112 


257 


0.205 


lib 


9 


6.79 


229 


0.20 


257 


0.205 



Table 2: The numerical values for Table[TJ These are taken from [11]. fj, = ^ 9 /(2ttT) with T = 0.2GeV, 
= v3/(27r) and fi g = \/3/id. Here we should notice that we will treat fj, g as a temperature-depending parameter 
which is different from the consideration in our previous paper [13] and will be discussed below in this chapter. 
The numerical values in this Table are taken from our previous paper [13]. 

Here we would like to touch on the field theory dual to our model (|3.2j) with the geometry (|3.1[) . 
Although it is not an exact statement as our model is a bottom-up model and does not stand on the 
configuration of D-branes, the assumed dual field theory would be D = 1 + 3 SU(7V C ) gauge theory with 
U(l) baryon symmetry and SU(Nf)L x SU(Nf)u chiral flavor symmetry as global symmetries at large 
't Hooft coupling at finite temperature, where Nf = 2 and the large- N c limit is taken, and no charge is 
concerned for the global symmetries in the model. 
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We have shown in [12] that the quark number susceptibility blows up when temperature is around 
160 ~ 190 MeV, which is related to the chiral symmetry breaking/restoration. We have further examined 
the mass spectra for the vector and axial- vector mesons as the function of temperature in ref. [13], which 
has suggested that the critical temperature locates at around 200 MeV. 

The temperature and the entropy density in the dual field theory are given from Hawking temperature 
and Bekenstein-Hawking formula as 



T H = ±- and s= ^2) = ^r^ 2 , 
irz 4G 5 2tt 



{(^T) 2 +^} 3/2 , (3.5) 



where we have used the relation G5 = ir/(2N 2 ). The relation of gu{i) an d 9su(2) to the soft-wall model 
are given as [9] 

2 16t 2 , 2 12tt 2 

9u{i) = -fp and 9su(2) = ( 3 - 6 ) 



The equation for the dilaton can be obtained from the equation of motion for the scalar field as 



$'(z) = 



3o'(z) , {f(z)v'(z))' a\z) 



m 2 x v{z) + -v 3 (z) 



(3.7) 



a(z) f(z)v'(z) f(z)v'(z) 

It will be noticed that in our previous paper [13], we have fixed X(z) with considering a regularized term 
Vi(z)kif(z) as X(z) = ln/(z) +vq{z) X I2 (and in accordance with this extra term, the coupling 

parameter A has to be modified), where vq(z) in this equation corresponds to v(z) in this paper. The 
explicit form for v\{z) is referred to ref. [13]. As it can be seen that without considering the extra term 
v\(z) \nf(z), the dilaton $'(w) in the vicinity of the horizon starts with the order (u — 1) _1 as 



u — 1 



with the numerators 
$_ 1 (T, M ) = - 



.4 



A (6C(ji)4 - 4) - 2Bz 2 {C(ji)4 + 6) L"\ m ^ ^ " 
+Bzl[m 2 x (4^ 2 + 1) (C(fi)z% + 1) + 2 {C{^ + 6) } 



(3.8) 

1) (C( P i)zt + l)-&C(p)zt + A) 

(3.9) 



where u = z 2 / Zq and [i = fi g /(2nT) (Later they are defined at eq. (|4.1ip and (|4.12p ). and A, B and C(/Lt) 
are given in Tabled] Then it turns out that the contribution $_i(T, fJ>)/(u — 1) appears at the order 
(u — l)~ x in the equation of motion for the fluctuations at the vicinity of the horizon. Therefore, the 
contribution $(^t) appears at the order (m — 1)~ 2 in the equation of motion. This prevents the solutions 
of the fluctuations from taking the in-falling boundary condition, where the in-falling boundary condition 
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is determined at the order (1 — u)~ 2 in the equations of motion. For this reason, it has been found in [13] 
to be useful for introducing the extra term vi(z)hi f(z) in order to regularize the divergence and obtain 
the in-falling boundary condition healthily. 

However, for our present purpose in this paper, it turns out that the extra term vi(z) hi f(z) makes 
the hydrodynamic analysis very complicated. Thus we have to discard the technique adopted in [13] with 
introducing the extra term v\(z) ln/(z) and find an alternative way. A simple idea is to consider the 
energy scale fi (or n g ) as a temperature-depending parameter n(T), so that one can adjust the parameter 
\x to make the numerator $_ 1 (T, \x) vanishes at the order (1 — As a result, we obtain the \i which 

vanishes the numerator <E>_i(T, /z) at each temperatures as shown in Fig[TJ and finally we can take the 
in-falling boundary condition. 
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Figure 1: We plot the positive and real values of fi found as the function of temperature. They are obtained from 
vanishing the numerator <I>_i(T, /j.) at the order (u— l) -1 in eq. (|3.9p . As a consequence, one can healthily take the 
in-falling boundary condition. The left and right plots are different just in the scale of x-axis. The lower curve (red 
points) and upper curve (blue points) in each figure represent two branches of solution for a given temperature. 



In FigJTl it is noticed that there are two branches for the solution, which are plotted in each figure by 
the lower curve (red) and upper curve (blue) for a given temperature. Then the question is which branch 
is physically meaningful. Before answering the question, we may make a comment on the relation of /z in 
our previous paper [13] as shown in the rightmost column of tableful and in this paper given in FigJTJ 

For the result of the model lib, two branches start at T c = 0.2 GeV, and the value of \x at the point for 
the branch of the lower curve (red points) matches also mostly the same value of \i shown in the rightmost 
column of tableEJ Then, let us consider which branches we should take. 

First, we shall recall that \i is the factor appearing in the factor a 2 (z) in ca. (|3.ip as a 2 (z) = l/z 2 + fx 2 = 
l/z 2 + (2irTfi) 2 . Then, one can see that if fi remains finite at high temperature, as temperature increases, 
our bulk space-time becomes completely different from the Schwarzschild AdS§ black hole space-time. It 
means that if /i does not vanish as temperature goes up, the gauge/gravity correspondence becomes invalid. 
On the other hand, if fi vanishes at high temperature, as temperature increases, our bulk space-time goes 



back to the Schwarzschild AdS$ blackhole space-time. At that time, since our model (|3.2[) has the extra 
terms (dilaton, scalar field and gauge fields) other than the Einstein-Hilbert action, the dual field theory 
is not the Af = 4 supersymmetric gauge theory even if the background geometry is the Schwarzschild 
AdSs blackhole space-time. However, what the background geometry goes back to the Schwarzschild 
AdS§ blackhole space-time in the high temperature limit would be a consistent condition as a holographic 
AdS/QCD model. In this sense, eventually, the physically acceptable branch would be the one represented 
by the lower curve (red points), and the rest would be inconsistent. 

An interesting point is that the effect of (i does not disappear abruptly, but gradually disappears as 
temperature increases. The effect of fi is a specific one in our AdS/QCD model, and it would be interesting 
to study its phenomenon in the future. Further, it would also be interesting that the branch of the model 
lib starts at about T = 0.2 GeV, where the temperature T c = 0.2 GeV is roughly the critical temperature 
for the Hadron/Plasma transition. 

4 Preliminaries for Hydrodynamic Analysis 

Before making a hydrodynamic analysis, let us first clarify the indices used in this paper 



M, N = x , X!, x 2 , x 3 , z or 0, 1, 2, 3, 4 (4.1) 

fi, v = x , xi, x 2 , x 3 or 0, 1, 2, 3 (4.2) 

k = x\, X2 or 1, 2 (4-3) 

a,f3 = xi, X2, x 3 or 1, 2, 3 (4.4) 

a, b, c = 1 ,2 ,3 (indices for SU(2) algebra) (4.5) 



where xo is coordinate of time, x±, x 2 and x% are coordinates of space on the boundary and z is radial 
direction of the bulk. 

We arc going to examine the retarded Green function for the scalar mode of f/(l) baryon current, the 
vector mode of SU(2) flavor current and the shear viscosity from the tensor mode of the energy-momentum 
tensor through AdS/CFT correspondence. For this purpose, we will analyze the hydrodynamics in the 
bulk gravity of the fluctuations, 



9mn ~ 


-> GmN — 9M N + hjtfN, 


(4.6) 


g MN_ 


-+G MN =g MN -h MN , 


(4.7) 


A M - 


Bqm = + Am, 


(4.8) 


r>a 

B M ~ 


-> B 0M = + B Mi 


(4.9) 
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where the background of the gauge fields are vanishing. Without confusing, we simply use the same char- 
acter for the total gauge fields and the fluctuation of gauge fields. These fluctuations are in linear order 
in equations of motion. In this paper, wc will treat the dilaton as the static background. 



We show the classification of fluctuations in the SO(2) little group for x and y-directions as 

Scalar mode: h 00 , h 3, h 33 , h zz , h 0z , h 3z , h k k and A , A 3 , A z , 

Vector mode: /ioi, h 3 i, h Z i and Ai, (4-10) 

Tensor mode: hij — 5ijh k k /2. 

One may see the correspondence of the above classification with the one given in Sec|2] 

In AdS/CFT correspondence, the fluctuations of the bulk gauge field A^ acts as the source for the 
i?-charge current J M on the dual field theory. Since it is known that i?-chargc can be regarded as an analog 
of the U(l) baryon number in AdS/CFT correspondence, we will consider that A^ acts as the source 
for the global U(l) baryon charge current (For example, see ref. [21]). Further, Bl 41 and -Br.^ also 
act as the source for the global SU(Nf)L x SU(Nf)n chiral flavor currents where Nf = 2 in this paper. 
The gravitational perturbations in the bulk act as the source for the stress-energy tensor T^ v on the dual 
theory. 

In our following analysis, for convenience we will define the normalized radial coordinate 

and the normalized frequency, momentum and the factor 

m = — , q = — and u= (4.12) 

Finally, let us mention the gauge fixing condition in our present consideration. Like other papers (For 
example, ref. [22]), we choose the axial gauge condition and the Landau gauge as follows 

A z = 0, B a z = and h Mz = 0, (4.13) 
9^ = and d»B Ltfl = d»B Rtfi = 0, (4.14) 

where the Landau gauge is imposed in the four-dimensional space-time on the boundary. 

5 Analysis on the U(l) baryon current 

We shall first perform hydrodynamic analysis in the scalar mode shown in Table. (|4. 10[) of U(l) baryon 
current. To begin with, we can write the equation of motion for the fluctuation of Am in terms of the field 
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strength as 



d M (^^T)e-^F MN (u)) = 0, 



(5.1) 



where g(u) = det §mn{ u )- We will carry out analysis in the momentum space by performing the plane 
wave expansion as 

' dw dq 



(27T) 



(5.2) 



where we have rotated the direction of momentum to x%, and pay attention to the same notation used in 
the gauge fields before and after Fourier transformation. 

cjg 00 (u)A' (u) - qg 33 (u)A' 3 (u) = , (5.3) 
d u (e-^^g^)g oa (u)g^(u)A' (u)) = e"*H ^gju) g ™{u)g 33 {u) (ujqA 3 (u) + q 2 A (u)) , (5.4) 



d u [e-^^g^u)g 33 (u)g uu (u)A' 3 (u)J = e"*M ^g~(ujg 00 (u)g 33 (u) (ujqA a (u) + w 2 A 3 (u)) , (5.5) 
du (e-^^glu)g aa (u)g^(u)A' a (u)) = ^gju)g^{u) (g™u 2 A a (u)+g 33 q 2 A a (u)) , (5.6) 

where ' = d u . 

From cqs. (|5.3[) and (|5.4p . the equation of motion for A' t can be obtained as follows 



(I) 



zq \ 2 d 
du 



d u (e~* W ^gjujg 00 (u)g uu (u)A' Q («) ) 



e^ W \/ I 9Wj 00 («)s 33 («) 
For a technical reason, we will perform a rescaling as 

2 



5 00 (u) 



S 33 (u) 



ro 2 + q 2 ^(u) = 0. 



-* {u) V^glujg 00 (u)g uu (u)A> (u) = A' Q (u), 



(5.7) 



(5.8) 



where we attach the factor — 2/zq in the front to make the leading term of the rescaling as unit in the 
expansion around u = 0. Note that when taking /i = 0, it recovers the case of Schwarzschild AdS§ black 
hole geometry. With the above rescaling, the equation of motion is given by 



(?) 



zq\ 2 d 
du 



d u A' {u) 



QV + q 2 



e -^ u ) y/-g(u)g 00 {u)g uu (u) 



A' (u). 



(5.9) 



e"*(") y/^g{u)g 00 (u) <? 33 (u) 
To make the hydrodynamic analysis, let us formally write down the solution to the first order hydro- 
dynamics as follows 

=,*(«) 



A' (u) = C^\u) 



= (1 - uT ir °/ 2 (l + toff^u) + q 2 G<%) + • • • ) 



(5.10) 



where Cq°\u), F^°\u) and G^(u) are going to be fixed below. It is noticed that the in- falling boundary 
condition has been used to yield the factor (1 — 

u )-ira/2 m 

the solution. This is realized by taking the 
special solution of fj,, as shown in FigQ] with requiring the numerator in the contribution at the order of 
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(u — 1) 1 to be vanishing in the dilaton &(u). Then we can see readily that Cg ' can be treated as a 
constant, 

cf (u) = C< 0) , (5.11) 

and the rest coefficients are obtained by the following equations, 

4 Q f iC™+ 2(1 -u)d u F?\u) \ = 
8 u \(l-u)e-^ u )^/^{u)g 00 (u)g 33 (u) J 

where the former and later can be obtained from the order at Id and q 2 in cq. (|5.10|) . Let us solve the above 
equations by taking the following three steps: 

1. performing the integration on the whole equations to get rid of the derivative d u {- ■ ■ } or d u (- ■ ■ ), 

2. rewriting the equations in the form "d u Fi = • • • " or u d u G^' =•••", 

3. carrying out the integration over the whole equations to obtain f[ ^ and Gf\ 



we then obtain the coefficients as follows 




2 , 



(5.15) 



with 



(0) _ qz 2 (4 U) q + 4 U ^) 



Here Cg ^ has been fixed by using eq. (|5.4[) with the Dirichlct boundary condition that the boundary value 
of Aq and A3 are given by A^ and A^ . Note that if the full integration in step 1 or 3 is difficult, instead 
of it, one can first carry out the Tyler expansion around u = 0, and after that perform integration toward 
its low order terms. As we will use eventually the GKP-W relation, the resulting solutions are the ones 
around the boundary. Where F^\ F^ \ G^' and G^ are the integral constants. F^ and G^ appear 
from the integration in step 3, and they can be taken arbitrarily. On the other hand, F^ and G^ appear 
from the integration in step 1, and they are in general fixed such that the solutions do not diverge at u — 1. 
Where the subscripts U H" and "B" denote the relevant constants which are fixed from the horizon and 
the boundary, respectively. 
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To put it more concretely, after step 2, one can see that F^'(u) and G^'(u) behave as F^'(u) , Gi'(u) 
l/(u — 1) around u = 1, which generally leads to the solution with logarithmic divergence, such a solution 
is ill-defined. To obtain a physically meaningful solution, a simple way is to vanish these contributions by 
exploiting the integral constants F^ and , so that the numerators become zero. 

In the actual calculation, F^ can easily be fixed to be 

4 0) = - ° ; e *° . (5.17) 

2yiT7? 

while fixing G# is in general complicated as the integrating becomes difficult due to the dilaton. For this 
reason, we have to perform integration around u = 1 to fix the integral constant Gh ■ It will be shown that 
as long as choosing the integral constant Gh appropriately, wc can arrive at the needed solution which 
has no logarithmic divergence and becomes well-defined in the whole bulk. Let us conduct step 1 with 
expansion around u = 1 

duG?\u) 



du- 



C, 



(0) 




e -*(") v / '-g(u)g uu (u)g 00 (u) 
du\c +C 1 (u-1) + C 2 (u - l) 2 • Cya L) 3 • •••} 



= G^> +C u + C 1 ^—-u) +C 2 [—-u' + u) +C 3 [ — -u- i + -u' 

= G H+f2^ C n + Cot + C 1 (^-t)+cJj-f+^+cJ (5.18) 



n=0 



with 



<7 (0) / 1 



and t = u-l, (5.19) 



e -*( u )^-g(u)g tl "(u)c/ 00 (u) / 

where "(n)" means the number of derivative with regard to u. In the above equations, we have performed 
an expansion from the second to the third line and written its result in symbolic way. From the third to the 
forth line, we have explicitly written the integral constant Gjj . From the forth to the fifth line, we have 
changed the variable u to t = u — 1. It is seen that when we rewrite it in the form that ll d u G^ = ■ • •", 
there is no term with l/(u — 1) in the expansion beyond the constant term. Therefore, we can readily 
obtain the expression of Gg to be 



G W = _yLiL,7 Iq *o G / {fW}), (5.2O) 



where we have rewritten symbolically the concise expression for convenience in the following actual cal- 
culation. Here G A (/i, {$ (ti) }) is a function with {$(")} = {$W(w)|„=i, $ (2) (u)|„=i, $ (3) («)I«=i, ■ ■ " } > 
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when n = and {$(")} = {O}, one has Ga(0,{0}) = 1. In this paper, G^ will be treated symbolically. 



It can be shown that the solutions near the boundary behave as follows 

4(«) = ~\ {G%\ 2 + F { >) Iog(«) + C (0) e $(0) (l + G^q 2 + ») , (5.21) 
= ^ f 4 0) q + ^7) log(-) - cf e"*^ (l + GfgV + Ffto) . (5.22) 
The boundary action So is given from the on-shell action in the quadratic order of fluctuation as 



5 ° = 33 / 7^ {[ 1 + ^3^ 4)(«K(") - 1 + 



)^4( U )4» . (5.23) 

/ a ) 



g 2 z 2 J (2tt) 2 V 20 / V z o 

By adopting the prescription for obtaining the retarded Green function in AdS/CFT correspondence [18], 
we arrive at the following two-point retarded green function of U(l) baryon charge current in the scalar 
mode 



_ N 2 T 2 e^-^V V + I* 2 



G W 8 2i^Tu~G {^{^)})k 21 (5 - 24) 

_ N 2 T 2 e^^VWTTc 2 
(Rl - — S — n- . T, 77. (O./iOj 



r(R) 8 2i7rTw-G (^,{$( ri )})fc 2 
From these results, we can read out the diffusion constant (12.91) as 



D= G A (,,{^})_ (5 . 26) 
6 Analysis on SU(2) vector and axial- vector flavor currents 

Let us make hydrodynamic analysis for the vector mode shown in Table (|4.10[) of SU(2) vector current 
and axial-vector flavor current. To begin with, let us combine the SU(2) gauge fields into the vector field 
Vm and the axial- vector field Wm 

Vm(u) = ^(B l . m (u) + B r . m {u)) and W M {u) = \{B L>M {u) - B R>M {u)). (6.1) 

The equations of motion for the vector field V£ and the axial- vector field W™ are given as 

V : D VM (yqffie-*™ F v ' MN (u)) = 0, (6.2) 
AV : D WM (^/=^) e -*W F^ MN (u)) + g 2 u(2) er^^g^)v 2 {u)g MN {u)W a M {u) = 0, (6.3) 

where v(u) has been defined in eq. (|3.3|) . and "V" and "AV" mean the field strength consisting of the 
vector and the axial- vector fields, respectively. F VMN {u) and Dv,m{u) are dictated as Fy MN (u) = 
d M V§(u) - d N Vfc{u) + if abc V^(u)V^(u) (f abc is SU(2) structure constant), and D VM = d M {u)-iV M {u) 
{Fw,m n (u) and Dw,m(u) are similar). It is seen that ea. (|6.2[) can be obtained from eq. (|6.3|) by dropping 

14 



the term proportional to Qgu^y Thus our main task here will focus on solving eq. (|6.3p . 



Now let us write down the equation of motion (|6.3[) to the linear order, and its result is given by 

d u (e-^v^)g aa (u)g™(u)W*i(u)) = e~^g aa (u) (g 00 ^ 2 + g 33 (u)k 2 - g 2 sui2) v 2 (u))W^u). 

(6.4) 

which cannot simply be solved by the double integral dictated above ea. (|5.14[) due to the extra term 
proportional to Sgj/^^M m the potential part. Thus we may factorize it as follows 

WZ(u) -> p{u)W a a {u). (6.5) 

Then cq. (|6.4p becomes the following form 

d u (e-*Wy/^)g aa (u)g uu (u)WZ'(u)) + p{u) {fi(u) - e~*M ^g P g aa (u)(g°W + ,g 33 fc 2 )} = 

(6.6) 

with 



Q(u) ee d u (e-^^^)g aa {u)g^{u)p\u))-g% u(2) e-^ (6.7) 

which indicates that when p(u) satisfies il(u) = 0, it becomes available to obtain the solution W°(m) from 
the double integral. Nevertheless, to fully solve the equation fl(u) — is technically difficult. Practically, 
considering that the needed solution at the last stage is the one only in the vicinity of the horizon and the 
boundary, we shall obtain the solutions in the expansion around u = and u = 1. 

First, let us obtain the solution in the expansion around u = 0. For this purpose, expanding Q(u) given 
in ea. (|6.7p . we arrive at the following form 

(ao + 0(u))p"(u) + (fa + 0(u)) P '(u) + + 7l + p(u) = (6.8) 

with 

a Q = \, /3 = -^($»| u=0 -2/i 2 ), 7o = \g% u{2) A 2 , 11 = g 2 su(2) A 2 (Bz 2 + 3Ap 2 ). (6.9) 

As a result, we obtain the general solution p(u) as a linear combination of two confluent hypcrgcometric 
functions ±Fi and U with two integral constants. We set two integral constants in such a way that \Fi 
vanishes and the leading of p(u) in the expansion of u = becomes 1. Eventually, we obtain p(u) satisfying 
eq. ()6.9j) to the expanded order as 

, w ^_^„) r(l _ ? ),(_|,„,A„) ( , 10) 
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with A = v /(/3 ) 2 -4a 7i. Where the functions T(u) and U(u) denote the gamma function and the 
confluent hypergeometric function, respectively. 

Now we shall try to solve fl(u) = in the vicinity of u = 1. It turns out that a general solution is given 
by a linear coupling of Confluent Hypergeometric Function U and Laguerre Polynomials L n . However, it 
finds difficult to fix the integral constants analytically, such that the constants in the solution at u = 1 
becomes common with the ones at u = of cq. (|6.10[) due to the lack of information for the solution of 
p(u) in the interior domain of the bulk. Therefore, we will only treat the solution p(u) around u = 1 
symbolically. 

Having obtained p(u) satisfying Q(u) — as in eq. (|6.10[) . we can solve the equation of motion by using 
the double integral similar to the previous chapter. With the results of p given in Table [1] we can formally 
write down the solution of W° in the hydrodynamic expansion with regard to ro and q 2 as 

W a a = C& o) (u)(l - u)- tm/2 (l + K>F[ a \u) + q 2 G[ a \u) + ••■). (6.11) 

With the same way stated above eq. (|5.14[) . the coefficients can be written as 

C ( a a) {u) = (6.12) 

and 

(op(a) $(o) ;r<( a ) \ 

?(o)=» -^rr+pfr 2 ). (fU3) 

(bg(u) - 1) + ZU £* Q) \u + 0{u% (6.14) 

where W< 0) denotes the boundary value of W£, and F^\ G^ and f£\ G% ] mean integral constants 
which are fixed at the boundary and the horizon as the same as the integral constants appeared in eq. (|5.14[) 
and (|5.15p . To be explicit, they are fixed as 

= liC™ v^?TIp 2 (l)e-*W, (6.15) 
G [ h ] =-\c^^Tie-^) Gw {p,{^}), (6.16) 

where G w (p,{¥ n ^}) is a function as the same as ea. (|5.20j) and Gw(0, {0}) = 1. From the above analysis, 
we can write down the solutions in the vicinity of u — with taking the factor p 

4p(0)e*(°) (f^u + G [ H ] q 2 ) + p 2 {0) (2q 2 Iog(«) - 2q 2 - ilu) 



2p{Q)[F < g ) u + G^q* + l 
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It is noticed that V£ can simply be read from the solution of W£ by taking g 2 su ^ 2 ) = 0- 

Finally, we come to consider the two-point retarded Green function in the dual field theory through 
GKP-W relation. For this purpose, the boundary action in quadratic order at u = is needed. One can 
see that it is given as 

which enables us to obtain the two-point retarded Green function in vector mode as 
lN 2 T 2 f 4p(0)e*(°) (F^uo + G^qA + p(0) 2 (2 g 2 log(«) - 2ry 2 - iu>) 




The two-point retarded Green function Gy a R (u) for V£(u) can simply be resulted from the G^v r( u ) given 
in above with vanishing g 2 SU (2)- 

7 Analysis on the tensor mode 

We are now going to carry out the analysis for the tensor mode of the fluctuation shown in Table (|4.10[) in 
the dual gravity, which allows us to calculate the shear viscosity (|2.7[) and its ratio to the entropy density 
in the dual field theory. For this purpose, we will start with the equation of motion for the bulk gravity, 
which can be generally given as 

Rmn — -GmnR + &Gmn = k 2 Tmn, (7.1) 

where k 2 = 8irG$ (G5 is given below cq. (|3.5[) ). A = — 6/Z 2 (I denotes the curvature radius of the AdS space 
and is taken to be 1) and Tmn means the energy-momentum tensor. From eq. (|3.ip . it is given as 



Tmn = g MN C + g PQ { -^—FpmFqn + Tr(F L ,pmFl,qn + Fr^pmFr^qn ) > — 2TyDmXDnX 

, 9u(i) 9su(2) v ' J 



(7.2) 

where C means the Lagranigian (|3.2[) except for the Einstcin-Hilbcrt part. It is noticed that when substi- 
tuting the background (|3.1j) into the above equation of motion for gravity, there is a deviation proportional 
to fig, which indicates that the background (|3.1[) is not an exact solution of the equation of motion for 
gravity as the back-reacted geometry is not considered here. Nevertheless, the effect from the back-reacted 
geometry was found to be small [20], and also it turns out that the deviated part appears only in the 
diagonal part of the equation of motion, while in our present treatment it involves only the off-diagonal 
part of the equation of motion, (x, y)-component, so the deviated part will be discarded in this paper. 
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To evaluate the viscosity, we will consider (x, y)-component in ea. (|7.ip . Then the energy- momentum 
tensor at the linear order becomes as T xy = Tr (g wu \d u X\ 2 — m 2 x \X\ 2 — j\X\ A } h xy , and the equation of 



motion for h xy can be obtained as 



= K y (u) + g j%Ky{u) + 92{u)h xy {u) (7.3) 



with 



_ y/uf(u) 

5l(u) = ^T' 

ff2(w) " 53(u) + 4^F ( " 2 - fc2/(u)) 

93 («) E 



2ua(u)f(u) 



lQuf(u)a"{u) + 8a'(u)(2uf'(u) + f(u)) + a{u)(s2uf"{u) + f'(u)) 

(7.4) 



-64a(u) 3 + K y^ ) T xy (u) 



Auf{u) 

Here (73(14) represents the constant part toward to and k in the potential term gi (u). It is seen that 53(14) is 
the obstacle to utilize the double integral dictated above eq. (|5.14p . which is similar to the case of eq. (|6.4p . 
Thus we may adopt the same analysis as the one for cq. (|6.4[) . namely exploiting the factorization, 

hxy(u) = x{u)<p{u), (7.5) 

and vanishing the extra term (73(14). 

Under the factorization (|7.5[) . we can rewrite eq. (|7.3l) as 

(ffi(u)xW(«))' + 9Mx 2 {u) U{u) + ( W 2 - k 2 f{u))^(u) = (7.6) 

with 

*(«) eee X "(u) + 9 -^\x'{u) + g 3 (u) x (u). (7.7) 

Then xi u ) is obtained by requiring ty(u) to satisfy ty(u) = 0. However, it can be shown that there is 
difficult to solve the equation ^(u) = fully. Similarly to the previous chapter, as the needed solution 
eventually is the one in the vicinity of the horizon and the boundary, we may consider the solution of x( u ) 
only in the vicinities at u = and u = 1. Nevertheless, from the reason mentioned below, one cannot 
obtain analytically the solution in the expansion around u = 1. 



Then we will try to obtain the solution of x( u ) m the vicinity of u = 0. To this purpose, expanding 
the constant part 'J(u) around u = 0, we will consider the following equation 



-^ + i( 40 ^ 2 + K2r *>)L ) 



+i{ - 8 + « 2 WT> y (u)\ u=0 + T>: y {u)\ u=Q ) } 



X (u) = 0. (7i 



18 



To obtain the solution xi u ) nl the vicinity of u = 0, we may assume the following form as the solution of 

\(u) = +Xo +XiW + XiUlog(u) with x-i = l/z%, (7-9) 
u 

where xo> Xi an d XL are the coefficients to be determined further. While we have fixed X-l = l/ z o- It 

can be seen that if we assume the form of solution as in eq. (|7.9p . then x-i can in general be arbitrary (and 

each coefficients except for x-i can be fixed so as to satisfy eq. (|7.8p ). Thus one should fix x-i at the first 

place as shown in ref. [19]. By doing so, one can see that xi u ) plays the same role as the background metric 

1 1 4/i 2 

9aa at the vicinity of u = 0, where g aa (u) in the vicinity of u — behaves as g a a{u) = — — I g - +0(u). 

z Q u z 

As a consequence, we can put the source of </>(it) as 

h^ v x {u) with h^ v x {u) the boundary value. 
By substituting x( u ) with the form in eq. (|7.9p and solving the equations in each order of it, we can fix 
the coefficients as 

Xo = ^ (4V + K 2 T' xy {u)\ u=0 ) , (7.10) 



Xi 



2 

X-i 



24 (36m- + ^»L ) 10S(U) (36/i2 + K2T - (U) I-) { 1728 ^ + 92 ^>)Lo 
+«2 (WL =0 + «a (^(«)L =0 ) 2 ) } - 2^2 {6T^( U )| u=0 (268/ + - 2 

+3T^( U )| u=0 (54 M 2 + K2T' xy {u)\ u=0 ) + 27^(0)} - 16128^ 



(7.11) 



XL = {l728 M 4 + 92 K2M 2 r; y (0) + « 2 (tT»)L =0 + « 2 (T^(n)| u=Q ) 2 ) } , (7.12) 

where we have obtained the solution of x( u ) to the order of u, which will become necessary in the analysis 
in what follows. 

Now we come to discuss the solution x( u ) m the vicinity of u = 1. From the expansion of the coefficients 
of cq. (|7.8[) around u = 1, one can see that a general solution in the vicinity of it = 1 is given as a linear 
coupling of Confluent Hypcrgcometric Function U and Laguerrc Polynomials L n . However, similar to the 
case of p(u) in the previous chapter, it is difficult to fix the integral constants analytically, so that the 
constants are common between the ones at u = 1 and u = (eq.([7jn]) with the coefficients (JTTTUJ) ) , due to 
lack of the information of the solution of x( u ) m the interior domain of the bulk. Thus x( u ) around u = 1 
will be treated symbolically in the present consideration. 



In the same way with the previous sections, we can obtain the hydrodynamic solution for the off- 
diagonal component of the fluctuation h xy (u) in cq. (|7.6p with xi u ) determined above. Consequently, it 
can be written as 

h xy {u) = X {u)4>{u) = C<*> X (u) (1 - M )-' iro / 2 (l + t»fi ( %) + q 2 G[*\u) + • • • ) (7.13) 
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with 



a 



(<t>) 



iTg^Tf^V 



rli+r -- 



4' r 



G { *\u) = G 



2x-\z \og{u) (xo - M 2 X-i) + X-ifo (xo ~ M 2 X-i) 
2(x-i) 2 ^o 



G 



2(X-i) 2 4 



(7.14) 
(7.15) 

(7.16) 



Where * has been mentioned above eq. 10[) . Fjf , Fjf and G}J , Gg represent the integral constants 
fixed at boundary and horizon, which are similar to the integral constants appeared in cq. (|5.14[) and (|5.15|) . 
More explicitly, they are fixed as 



1 



v/VTi 
1 



u H - - 



(7.17) 
(7.18) 



V V + 1 

Here we have used x—i = l/ z o- It enables us to write down the solution of h xy (u) in the vicinity of u = 

as 



h X y(u) 



h(0)x U0)x 

^^ i + ^(Xo-q 2 X-i) + A ^ 



+2zfa 2 X - 1 (xo - M 2 X-i) Iog(u)l u + G(u 2 ) 



i 4^ + (? 2 {G^+x-i4(xo-M 2 X- 



(7.19) 



which shows that the expression diverges at u = 0. However, such a divergence arises from the contribution 
of x( u ) given in ea. ([7.10[) . and the solution of (/>{u) itself has been obtained healthily as it can be seen 
from eq. (|7.13[) . Such a situation is the same as the one discussed in ref. [19] and also other papers based 
on ref. [19]. 



Now that we have obtained the solution, let us evaluate the viscosity through Kubo formula. The 
boundary action So is given by 

5*0 = Son-shell + SgH + S c t, (7.20) 

with 

Sgh = j^- [ d^xyf^K, (7.21) 

07TG5 J 

where Sqh and S c t denote the Gibbons-Hawking term and the counter term respectively, and gpO) is the 
four-dimensional induced metric, and K and RS A ' mean the extrinsic curvature and the curvature on the 
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boundary toward our geometry (|3.ip . respectively. They are found to be 



Son-shcn = ttt- I CTe-* (u) 4^|4a'(«)/(u)/i*>) 2 - 3a(u)f(u)h xy (u)h xy '(u)} , (7.23) 

2k z zo J (27r) z a{v,y I J u=o 

Sgh = Um / ^e-^^{h xy (uff'(u)+4f(u)h xy (u)h xy '(u)}, (7.24) 
2k z zo u-¥0 J (27r) z a(w) I. J 

1 /• rfwdg ^_ $(tt) 12a(M) 2 /( M ) + fc 2 /(M) -cj 2 2 

5 ct = — hm I -^zv2 e ; 7-77 — f KyW- I 7 - 25 ) 



8 K 2 u^oy (2tt) 2 ' a( W ) 2 y/M * y 

The retarded Green function is obtained by using the GKP-W relation 

N 2 T 2 



GM^u, k) = -^-{c^fc 2 + 27TT (ii#> w - x 2 -x) }• (7.26) 



which allows us to evaluate the shear viscosity via Kubo formula 



lxa(Giy > X y(ui,k = Q)) 

7] = - lim ^ '- = -F H N 2 T 3 . (7.27) 

c^o uj 8 



with this result, it is then not difficult to obtain the ratio 

V 1 



(7.28) 



s 47iVl + V(l + 16tt 2 M 2 T 2 ) 3 / 2 
where the entropy density is given in ea. (|3.5p . 

The ratio rj/s as the function of temperature is shown in Figj2l where the left plots and right plots is 
just in the scale of s-axis. The two curves (upper red points and lower blue points) represent the results 
obtained by using the ijl in FigfTJ and the each plots in FigJ5]and Fig [1] correspond each other. The dashed 
horizontal line represents the KSS bound, 1/(Att) [23]. It can be seen that there are two branches with 
regard to the solutions of fi. As it has been discussed in the ChapO the physically acceptable branch is 
given by the lower curves (red points) in the plots. 

8 Conclusions and Remarks 

In this paper, we have extended an IR-deformed AdS/QCD model presented in ref. [11] to finite tempera- 
ture system by another way different from the one proposed in ref. [13], worked out the hydrodynamics and 
computed transport coefficients. In this extension, we have found several branches, and we have chosen 
the branch in which our model can be consistent as a holographic model. 

The specific things in our model are the parameters /i g , gsu(2)j A the scalar field X and dilaton $. 
Among these, the ones entering to our analysis in this paper has been [i g thorough /1 = fj, g /(27rT) and 
dilaton. We have treated the effect of these in symbolic way. The rest things will enter if one considers 
the background dual to the system at the finite chemical potential. However, the analysis with the finite 
chemical potential might be too complected to perform the hydrodynamic analysis. 
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Figure 2: We plot r\/s given in eq. (|7.28[ ) using given in Fig[T] The lefts and right plots are different just in the 
scale of x-axis, The red and blue points represent the results obtained by using the in Fig[T] These plots and the 
plots in Fig [T] correspond each other. The dashed horizontal line represents the KSS bound, l/(4.7r) [23]. Among 
these plots, the result represented by the red points corresponds to the results concerning the /j, we have chosen as 
physically acceptable one. 

A further direction of this study will be the holographic description of Bjorken flow [24]. Bjorken 
Flow is an effective model having been invented to describe evolution of quark-gluon plasma produced in 
high energy collision experiment. This has the property of the ideal fluid, and diffuses with the form of 
distribution in the boost invariant way from a point in Minkowski space-time where a collision happens, 
it is considered that the quark-gluon plasma being produced and diffusing in RHIC or LHC, etc can be 
described by using Bjorken Flow. However, it is known that understanding of Bjorken Flow in the frame 
work of the field theory is hard due to its strong coupling and real-time evolution. In such a situation, 
the gauge/gravity correspondence would be an effective measure. For the relevant studies, we would 
like to refer reader to the reference in the review papers [16]. In the analysis of these, we can see that 
hydrodynamic analysis is necessity. Therefore, using the hydrodynamic analysis formulated in this paper, 
we are going to try to the Bjorken flow in our model. 
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